Transport through a quantum dot with excitonic dot-lead coupling 
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We study the effect of a dot-lead interaction on transport through a quantum dot hybridized to 
two semi-infinite Luttinger-liquid leads. A bosonization approach is applied to treat the interaction 
between charge fluctuations on the dot and the dynamically generated image charge in the leads. 
The nonequilibrium distribution function of the dot and the tunneling current are computed within 
a master-equation approach. The presence of the excitonic dot-lead coupling is found to enhance 
transport in the vicinity of the Coulomb-blockade threshold. This behavior is in contrast to the 
usual power-law suppression of electronic tunneling which is found if this interaction is ignored. 
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I. INTRODUCTION 

Understanding nonequilibrium phenomena in trans- 
port through quantum dots and single-molecule devices 
is of fundamental interest and is a major challenge in 
the field of nanoscience. Much attention has been paid 
to the interplay of the dot-lead hybridization and on-dot 
interactions involving vibrational degrees of freedom or 
the local repulsions which lead to magnetic moments and 
to Kondo physics. However, the effect of an 'excitonic' 
interaction between the charge on the dot and the charge 
on the leads has been studied little, even though it clearly 
should be present on physical grounds. In this paper we 
address this physics in the context of a quantum dot con- 
nected to two one-dimensional leads. The model we con- 
sider is inspired by recent experimental investigations of 
the conductance of a molecule placed in the gap created 
by breaking a carbon nanotubeji^ although theory and 
experiment are not yet at the point where a quantitative 
comparison can be made. 

One may view a quantum dot with one-dimensional 
leads as an impurity embedded in a Luttinger liquid. 
Fabrizio et al£ and Maurey and Giamarchi^ have studied 
the case of an impurity described as a short-ranged po- 
tential scatterer without dynamical charge fluctuations. 
Lerner et al^ have taken into account finite dot-lead hy- 
bridization but have not considered a Coulombic dot-lead 
coupling. Goldstein et al& have investigated the effect of 
a Coulombic dot-lead interaction on the dynamics of the 
population of a quantum dot using density matrix renor- 
malization group and classical Monte Carlo simulations. 

In a previous paper we studied the effect of a cou- 
pling between dot and lead charge densities on the re- 
laxational dynamics of a quantum dot side-hybridized to 
a Luttinger-liquid leadii We found that this 'excitonic' 
coupling has important consequences for electronic cor- 
relations and may enhance the tunneling of electrons in 
the regime of weak hybridization. The present paper ex- 
tends our analysis to the case of a quantum dot placed 
between two leads such that it cuts the Luttinger liquid 
into two semi-infinite quantum wires. This situation dif- 



fers in several respects from that considered in our pre- 
vious work. Most fundamentally, intrinsically nonequi- 
librium behavior driven by a current flow across the dot 
is possible. Further, because we deal with semi-infinite 
systems, boundary exponents appear instead of bulk ex- 
ponents. Finally, the presence of two leads means that 
tunneling into one lead is modified by orthogonality ef- 
fects arising from the excitonic coupling to the other lead. 
These differences turn out to produce significant changes 
in the results. 

The rest of the paper has the following structure. In 
Sec. [XT] we introduce the model. In Sec. IIIII we explain 
the bosonization scheme which we use. In essence the 
idea is to write the model in a basis corresponding to a 
translation-invariant system and to note that the break- 
ing of translational invariance corresponds to imposing 
boundary conditions on the lead wavefunction which may 
be satisfied by choosing particular linear combinations of 
wavefunctions in a way inspired by the "method of im- 
ages" in electrostatics^— Section |IV] solves the resulting 
the model by canonical transformations a la Schotte and 
Schotte^ Results for the tunneling rates and current- 
voltage characteristics, calculated to leading nontrivial 
order in the dot-lead hybridization, are presented in 
Sec. [V]and Sec. ED Section Ivm presents a summary and 
conclusions, and indicates possible directions for future 
research. 



II. MODEL 

Our system consists of leads, a quantum dot, and hy- 
bridization and excitonic coupling terms. It is described 
by a Hamiltonian of the form 



H = H\ 



lead 



dot 



hyb- 



(1) 



For simplicity, we take the quantum dot to have one 
(spin-degenerate) level and a repulsive interaction. How- 
ever, we emphasize that the method is general and would 
work in more complicated situations. The quantum-dot 
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FIG. 1: Schematic of a one-dimensional lead-dot-lead sys- 
tem showing two leads, possibly maintained at a nonvanish- 
ing voltage difference, the quantum dot (filled circle, green 
online) and a gate which screens the long-ranged part of the 
Coulomb interaction. 



Hamiltonian fldot is thus 

fldot = £d n d + -jTid(nd - 1). (2) 

Here U is the local Coulomb repulsion, n d = d\d a is 
the operator giving the total number of electrons on the 
dot and d\ creates an electron with energy e d and spin 
(7 on the dot. 

We model the lead part of the Hamiltonian as two 
disconnected half-chains. Motivated by the possibility 
of nanotube leads we allow the two leads denoted by 
a = L,R to have an orbital (channel) degeneracy la- 
beled by 1 < a < N as well as a spin degeneracy. The 
lead Hamiltonian is then 

fllcad = £gkC niakrT C aa ka + -Hint- (3) 

ctakfj 

Here c^ aakrT creates an electron in a state of energy e ak 
in orbital a of lead a, and flint labels the interactions 
among lead electrons, which will be discussed in more 
detail below. 

The situation we consider is sketched in Fig. [T] The 
coordinate system is defined such that the quantum dot 
is at x = and lead a terminates at x = x a with xl < 
and xr > 0. We assume that in the absence of the dot 
the leads are decoupled so that the electronic boundary 
condition is 

^aao-(Xa) = 0, (4) 

where ip aa <j(x) creates a lead electron at position x. 

We restrict ourselves to the case of short-ranged in- 
teractions between the conduction electrons, assuming 
in particular that the interaction between electrons in 
different leads is negligible. The effect of interactions be- 
tween leads is an important open question. Interactions 
become short-ranged in the presence of a metallic gate, 
which screens the long-ranged part of the Coulomb force. 
We shall need an explicit form for the charge- channel in- 
teraction in the leads, which we take to be 

H ^ = \ E V « [ a dx[p*°\x)] 2 + ..., (5) 

a=L,R, J 



where J indicates an integral over the region of lead a 
and 

pT(x) = ^2p aa<T (x) (6) 

a a 

with 

Paaa(x) = ^ C aa{k+q)<r Colaka CO 
qk 

denotes the total charge density (summed over all chan- 
nels and spins). We shall not need to specify other po- 
tential contributions to flint- 

We include two forms of dot-lead couplings. The first 
is the standard hybridization 

flmix = E / dx E t 7 "^) + h - c -] > ( 8 ) 

a aa 

with the tunneling amplitudes T a {x) strongly peaked at 
x in the vicinity of the lead edges. The second term is a 
(screened) Coulombic dot-lead coupling, 

flexe = n d [ dxW a (x)p^ (x) , (9) 

where the dot-lead interaction W a {x) is also peaked for 
x near the lead edge. 

III. BOSONIZATION SCHEME 

A. General formalism 

Because the leads are one-dimensional and the disper- 
sion may be linearized near the Fermi level, we expect 
to be able to represent the low-energy degrees of free- 
dom in terms of bosons representing propagating charge 
and spin fluctuations^^ The semi-infinite nature of 
the leads presents a technical difficulty in that the stan- 
dard bosonization formulas are derived for translation- 
invariant systems, whereas in the model of interest trans- 
lation invariance is strongly broken. To deal with this 
situation we use ideas introduced by Kane and Fisher^ 
Eggert and Affleck^ and Fabrizio and Gogolin^ For each 
lead we introduce a translation-invariant reference sys- 
tem [Eq. ([3]) with the range of x extended from — oo 
to +oo], which we bosonize as usual. We then im- 
pose a boundary condition which selects from the refer- 
ence problem only those states which fulfill the physical 
boundary condition. 

B. Bosonization of reference problem 

We now recall the bosonization of the translation- 
invariant reference problem, in order to define notation. 
We combine spin and orbital quantum numbers into a 
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superindcx ft = 1, 



p^a{x) describing right (A 



, M and introduce the operators 
) and left (A = — ) moving 
particle-hole pairs in orbital /3 of the reference system 
corresponding to lead a. From these we define boson 
operators 0° j/9 , 6°^ via 



ve%(x) = - [pip{x) - p -p{x) 

that obey the volume commutation relation 
[4>%{x),V0 a al p,{x')\ = l -5 aa ,8 f3 f 3 ,8{x-x' 



(10) 
(11) 

(12) 



and are related to the total particle density in state (3 by 

p a p(x) = V<j)° ap (x). 
Expressing the lead electron creation operator 

A=± 

in terms of these bosons fields gives 



(13) 
(14) 



^ ( x ) = ^e^KW-Wi 



(15) 



Here the small positive infinitesimal factor r\ arises from 
the correct normal ordering of the operators^. The Klein 
factor U a carries the Fermi statistics and a time depen- 
dence related to the lead chemical potential p, a , 



U Xa (t) 



J Vat 



U Xa (0). 



(16) 



Under reasonable conditions (no Umklapp scattering, 
standard interactions) the low-energy part of the Hamil- 
tonian for the reference translation-invariant problem 
may be diagonalized in terms of boson fields </> 

ab j @ab 

which are linearly related to the original operators 
(f)° a p ,0^p and obey the same commutation relations, 

Eq. (fT2|) . 13 i 14 It is convenient to rearrange the 8 a i,, <^ a & 
into chiral fields 4>^ b conventionally defined by 



0± b (x) = 8 ab (x)±4> ab (x), (17) 
in terms of which the reference Hamiltonian becomes 

TTVab 

xb 



dx[[V<t>+ b (x)¥ + [V<t>- b {x)] 2 ) (18) 



with the velocity parameters v ab and coefficients of the 
linear transformation relating the <jy^ b to the <^g, O^p 
determined by the bare velocities and interactions of the 
lead cigcnstates. Note that [ffrab'&a'b'] = ^ while 

&(^).V^, b ,(y)l = ±-5 aa ,5 w 5{x-y), (19) 

J 7T 

which among other things implies that (jy^ b {x,t) are re- 
spectively functions of x ± vpt only^ 



In the general case the transformations arc compli- 
cated, but in the physically most relevant case where all 
channels in a given lead have the same bare velocity and 
the interactions conserve the total lead density, two sim- 
plifications occur. First, the transformation does not mix 
the 8 and 6 variables so that 



b 



6 > 



(20) 
(21) 



Second, one channel (which we take to be b = 1) describes 
the total charge density in lead a so that V(f> a .b=i is 
linearly proportional to p^ [Eq. ([6|)]. Specifically, 



,6=1 



.6=1 



J a.b=l 



KM 



and, conversely, 



9° 

y a/3 



1 IKr 

2 
1 
2 



(^,6=1 



IJ a,b=l 



(22) 
(23) 

(24) 
(25) 
± 



6>2 



with the ellipses representing the other terms </>^ 

(f>~ b>2 needed to make up the full operator. 

The 'Luttingcr parameter' K c and charge-channel ve- 
locity v a ,b=x = v c are related to the bare charge-channel 
interaction V£, Eq. (|5|), and the bare Fermi velocity vf 
by 



K r = 



1 



1 



MV„ C 



MVl 

TTVf 



(26) 



Transcribing Eq. (TT5"j) into the new basis yields 
il>X»p(x) = U Xa e lXkFX ^ x %{x) 

^ ( ^ VT?7 , _ 



x e 



(27) 



with i/'ao^ M an exponential of a combination of the cjr^ b 
with b = 2,...,M. 



C. Implementation of boundary conditions 

Equation (|14p. in combination with the constraint, 
Eq. (j4]), on the lead wave function implies 

U+ a e ikpXa e™^*^ + C/_ a e~ lkFX " e-^V*-*) = 0. 

(28) 
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for all P and a. Equation (|28|) is seen to be satisfied if 
we impose the two conditions 

(i) cos [7T0°g(x a ,i) + k F x a ] = 0, (ii) U + . a = U- >a . 

(29) 

Condition (ii) is the statement that perfect reflection 
at the channel edge means that the Klein factor, which 
carries the position dependence only via the chemical po- 
tential, is the same for left and right movers in each lead 
a so that henceforth we drop the subscript A in U\ a . 

Condition (i), in combination with Eqs. (f2"4"f — (|23|) 
says that at the lead boundaries x = x a the difference 
between the + and — chiral fields must be time indepen- 
dent. Shifting the origin of the coordinates for lead a to 
x a we find 

&(0,i)-<fe fe (0,*) = A*6 (30) 
with the constants D a b such that [using Eq. ([20]) ] 

Because the chiral fields ^ b (x,t) are respectively func- 



tions of x ± v F t only we can extend Eq. (l30l) to all space 
and time as 

0+ 6 (s,t)-^ 6 (-X,t)=X>a6. (32) 

Inverting Eq. (|3Tj) . noting that the left-hand side is inde- 
pendent of /? and using the fact that the charge channel 
is the same in original and eigenstate variables we find 

Dal = ^U- 2 J^\ D atb ^=0. (33) 

To compute correlation functions for operators cor- 
responding to the semi-infinite Luttingcr liquid one 
rewrites a general expression in terms of the chiral oper- 
ators <fy^ b , uses the boundary condition Eq. ([30]) to elim- 
inate, say, (j)~ b (x,t) in terms of <p^ b (—x,t) and D a t, and 
then computes the <\>\ b correlation functions in the usual 
way from the <f>* b part of Eq. (fT5|) . For example, the 
fermion operator becomes 



IV. ELIMINATION OF DOT-LEAD COUPLING 

The dot-lead interaction H cxc in Eq. ([9]) can be elim- 
inated by a canonical transformation as first noted by 
Schotte and Schotte.-^ Using Eqs. ©, ^T3j) and ([32]) gives 



\fKcM 



/ dxW a (x) 



v< fe=1 (*)-v>+ (-*) 



n d - 



(35) 



We now formally extend the integral over the full 

translational-invariant reference system by defining 

W a (— x) = W a {x) (recall we have shifted the origin of 
coordinates in lead a to x a ) obtaining 



#cxc = ^P^]T / dxW a (x)W^Ax)n d . (36) 



Examination of Eqs. (|18|) and (f35|) then shows that the 
dot-lead coupling may be removed by a shift V</>+ 1 (x) — > 
^4>ti( x ) ~ n dZ a (x) with 



Z a (x) 



VK^IW a (x) 



(37) 



After this shift the Hamiltonian retains its original 
form, except that H cxc has been eliminated and the dot 
parameters £d and U are shifted to Sd — > id = £d — A 
and U — > U = U — 2A with the polaron shift A defined 
by 



(38) 



Observe that the magnitude of the renormalizations de- 
pends on the magnitude of W a , which in turn may rea- 
sonably be expected to vary with the distance between 
the dot and lead a. 

The commutation relation Eq. (|I9[) shows that the mo- 
mentum conjugate to V0* h=1 is — 7r<^ fc=1 /2 so that the 

shift is effected by the canonical transformation O — > 
e iS 0e -iS with 

S = n d J2^ / dxZ a (x)4>+ b=1 (x). (39) 

Under this transformation the hybridization term Eq. ([8]) 
becomes 



e iS H ■ e~ iS 

° 1J mix t - 



£ J dxT a (x)4e 1 ^ 



B a (x) 



x e 



^T, a > fdx'Z a ,{x')4>+, b=1 (x') 



i>\ af s{x). (40) 
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Here the first factor in the second line comes from trans- 
forming the operator d\ and the second from transform- 
ing the lead operator. The constant B a (x) is given by 

B a (x) = A-y^ J dx' Bgn(aj') [Z a (x' + x) - Z a (x' - x)} 
i r 

dx' sgn(x') [Z a (x' + x) + Z a (x' - x)] . 

(41) 



We shall be interested in \x\ small compared to the range 
over which Z a is nonvanishing, in which case because W a 
is defined as an even function B a will be negligible. 



EXCITONIC INTERACTION AND 
DOT-LEAD HYBRIDIZATION 



The transformations introduced in Sec. II VI remove the 
explicit excitonic dot-lead coupling from the Hamilto- 
nian, at the expense of adding operator content to the 
hybridization. The resulting model can (to our knowl- 
edge) not be exactly solved, but it can be studied by 
perturbation theory in the dot-lead hybridization. We 
present here an investigation of the lowest nontrivial or- 
der in perturbation theory, which reveals the essential 
physics introduced by the excitonic dot-lead interaction. 
The methods introduced in Refs. (l6lfl7| can be used to 
extend the calculation numerically to all orders. 

In a pcrturbativc analysis of Eq. (|4T)|) the crucial quan- 
tities are expectation values of the form 



F a (t)= Y, (£JU(*)&'a/j(0) 



(42) 



X,X'=± 



with 



&a*» = U a / dxT a (x)e^ B ^ 



x e 



} E„/ /d*X,(*')<^,, 6=1 (*') ^ (a,). (43) 



Equation (|42|) is evaluated by using Eq. (|27|) to express 
the fermion operator in terms of chiral boson fields and 
then using standard results of bosonization. To avoid 
inessential complications we specialize to the physically 
most relevant case of short-ranged dot-lead interactions. 
In this case we may neglect the B a (x) term in Eq. (|43p. 
We will also assume that (as in the nanotube case) all of 
the non-charge channels are characterized by very weak 
interactions, so that we may assume that in all channels 
except the charge channel we have free-fermion correla- 
tions. In the long-time limit we then obtain (including 
for later convenience a nonvanishing temperature T) 



F a (t) 



"F (t)e* 

T n \Vr. ' 



(44) 



with Fo(t) = TrrjT/[vF sinh(7rTt)] the free-fermion corre- 
lator, 

— = 2tt / dxdx'T:{x)T a {x') 

T a J 

x 4 sin(fcF {x — x a )) sin(fcF(x' — x a )) (45) 
proportional to the bare tunneling rate, and 



(i-z a ) 2 



IV .A. 

x log( -^sinluVTt) 



(46) 



where we have defined a = L if a = R (a = R if a = L) 
and Z a = J dxZ a {x). Here A is the momentum cutoff 
of the Luttinger-liquid behavior; in the nanotube case it 
is in essence the tube diameter, which is the length scale 
on which the Coulomb interaction is cut off. 

The integral in Eq. (|4"5)) reflects the vanishing of the 
lead fermion operator precisely at the edge of the lead, 
x a , and may give rise to a strong dependence of electronic 
correlations on the position of the quantum dot (relative 
to the leads) but need not be evaluated precisely for our 
subsequent considerations. 

The main physical content of Eq. (|46|) is that the func- 
tion F a (t) is renormalizcd from the free-fermion behavior 
cx [sinh(TrTi)]- 1 to [sinh(7rTi)]" Y ° with 



Y a = 1 - — 
M 



1 - 



(i-z a ) 2 Z\ 

K r K r 



(47) 



Differences from Y a = 1 corresponds to changes from the 
free-fermion situation. Y a > 1 corresponds to a suppres- 
sion of tunneling and Y a < 1 to an enhancement. 

In the absence of the excitonic coupling (Z a = Z & = 0) 
we obtain Y a = 1 + (1 — K C )/K C M, the standard result 
for tunneling into the boundary of a Luttinger liquid with 
only density correlations! 15 ' 18 ' 19 The excitonic coupling 
has two effects. The (1 — Z a ) 2 /K c term leads to a weak- 
ening of the charge-channel renormalization and thus to 
an enhancement of tunneling, similar to that found in 
our previous workii However, the term Z\jK c , which 
did not occur in our previous work, leads to a strength- 
ening of the charge-channel renormalization and hence to 
a suppression of tunneling. This term arises because tun- 
neling into one lead changes the dot charge density. The 
excitonic interaction means that this change in the dot 
density causes an orthogonality catastrophe in the other 
lead, suppressing electronic tunneling. The total renor- 
malization is governed by a competition between the two 
effects, and is thus controlled by the relative sizes of Z a 
in the two leads as well as by K c . 

If the situation is symmetric and the interactions are 
of Coulombic origin with a reasonably large screening 
length, then as shown in Ref. one has Z a = 1 — K 2 
with K r < 1. In this case one obtains 



Y L = y R 



1 



1 + t 1 - Kc - 2i ^(l - K c)} ■ (48) 
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Especially in the strong-interaction case (K c <C 1), the 
value of Y a is changed relatively little from the Z a = 
value. For a nanotubc with a screening length some- 
what larger than the tube diameter, the considerations of 
Ref. imply K c on the order of 0.5-0.7. In this circum- 
stance, Y a is 1.1—1.3 in the absence of excitonic effects 
while in its presence Y a is 0.9-1.1. 

On the other hand, a very asymmetric situation (Zl — 
1, Z R ~ 0) would lead to Y L ~ 1 - 1/M < 1 and Y R ~ 
1 + (2/K c — 1)/M > 1, so that tunneling into one of the 
leads is suppressed and the other is enhanced; however, it 
is difficult to realize a small Y a even in the case of M = 1 
channel. 



(a) 




VI. CURRENT- VOLTAGE CHARACTERISTICS 



The assumption of weak coupling between the quan- 
tum dot and the leads suggests a master-equation ap- 
proach to study the nonequilibrium dynamics of the sys- 
tem. The quantum dot with a spin-degenerate level is 
described by the diagonal density matrix 



I (ek L ) 



Pd = Vc 



P 1 |l)(l|+P 2 |2)(2|, 



(49) 



where Vq, Pi, and Vi denote the occupation probabilities 
of the empty state |0), the singly-charged state |1), and 
the doubly-charged state |2). 

In the Markovian limit of sufficiently slow dot dynam- 
ics, inserting Eq. (|49[) into the von Neumann equation 
and expanding to lowest order in the hybridization yields 
the following set of master equations 

Vo = ^[PiKU-2W, (50) 

a 

A = ^[2W^i+2P2^ 1 -Pi(^o+^ 2 )], 



(51) 
(52) 



0.75 
0.5 
0.25 
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r=o.oi 



T=0.l 

r=o.3 



FIG. 2: Current-voltage characteristics of a broken Luttinger 
liquid (a) for Y = 0.9 (representative of the presence of the 
excitonic interaction in a symmetric junction with nanotube 
leads) and (b) for Y = 1.1 (representative of a dot with nan- 
otube leads in the absence of the excitonic interaction), and 
different temperatures. We show results for a symmetrically 
applied bias voltage eV = — with the gate voltage set 
to zero. We assume that the renormalized and bare dot en- 
ergies id and Ed are positive. Because our interest here is in 
the form of the threshold behavior we choose units such that 
in each panel the bias voltages V and thermal energies T are 
given in units of the corresponding onsite energy appropriate 
to that panel. The local Coulomb interaction is assumed to 
be twice as large as the onsite energy in both panels. 



with the tunneling rates 



U + p a ), respectively. Using Eq. (|44|) we find 



Tl^ = 2 Re / drF a (r) e~ lidT , 



(53) ft* oc 



ft?_+ 2 = 2Re/ drF a {T)e- l{Sd+u)T . (54) 
Jo 



The transition probabilities for the emission and ab- 
sorption of an electron are related by IZ^^^Ed — p a ) = 
Ki^ {-e d + p a ) and TZ^ 2 (e d + U - /j, a ) = fcgUi - 



- E ^'° /r, rri \ Y a -1 

e 2t / 2-kI 



2nT a \ v c A 



^ > 2 ^ 1 2irT 



(55) 



K?_,n oc 



e 



(2kT\ 

^ 2 2nr a \v c k) 



Y a -1 



2 ' 1 2nT 



(56) 



The rates in Eqs. ((55]) — obey the detailed-balance 
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(a) 



(b) 



FIG. 3: Differential conductance dl/dV (in units of its maxi- 
mum value) as a function of bias voltage V for different values 
of Y calculated for voltages in the neighborhood of the higher 
bias step shown in Fig. [2] at temperature T = id/ 30 and other 
parameters identical to those used in Fig. [2] 



V 



-4 







A V 





I. 



-2 2 

Vg 



4-2 2 
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FIG. 4: Color-scale plots of the current as a function of bias 
voltage V and gate voltage V a (a) for Y — 0.9 (representing 
the presence of the excitonic interaction) and (b) for Y = 
f.l (representing the absence of the excitonic interaction). 
Voltages are given in units of the onsite energy. We assume 
T = ed/100, U = 3e d , and A = e d /2 so that U = 2e d . 



condition, because each lead is individually in equilib- 
rium. 

The operator for the tunneling current through lead a 

is 

I a = ~ie^2 / dx [T a {x)d\ip aa<7 {x) - T*(x)ipl aa {x)d a ] . 

CUT 

(57) 

Solving Eqs. ([50 ]) — ([52 ]) in steady state (V t = 0) then 
gives the steady-state current 

IT \ — r 2 ^ (l? a T? & — 7? Q 7? a ") 
\i a / — e ^/vo->i' v i->o '^i->o'^o~>i; 

+ e (^ftf-i - fc§Uirc?_ 3 ) (58) 

with 7£„_y m = X) Q and 

i? = + 2110^112^1 + fto-n^i-*. (59) 

Thus we see explicitly that the steady-state current is 
proportional to the product of the hybridization to the 
left and right lead and changes when the chemical poten- 
tial in either lead lines up with one of the dot levels. In 
the asymmetric situation different terms come in to reso- 
nance at different bias voltages, leading to a diversity of 
resonance behaviors. 

In the low-temperature limit, the current-voltage char- 
acteristics obey a power law close to the Coulomb- 
blockade thresholds, eV = ±2e<j. Here the electronic 
correlation function is F a oc 1 /t Y " , which gives rise to 
tunneling rates of the form 

KSUi oc - ( ^M) Ya 1 e ( Pa - e d ) . (60) 

T a V VcA J 

Numerical results for the current-voltage characteris- 
tics for a symmetric junction (equal hybridization and 



excitonic couplings to left and right lead) are presented 
in Fig. [2 for exponent values representative of a dot cou- 
pled to carbon nanotube leads with (upper panel) and 
without (lower panel) excitonic couplings. The steps in 
the current are associated with the lead Fermi level com- 
ing to resonance with the n d = to n d = 1 energy dif- 
ference (lower bias feature) and the rid = 1 to rid = 2 
energy difference (higher bias feature). The symmetry 
of the situation implies that the steps are the same for 
positive and negative bias. We see that in the excitonic 
case (upper panel) the conductance is enhanced at the 
threshold, while in the absence of excitonic coupling the 
conductance is suppressed at the threshold. As is seen 
more clearly in Fig. [3] the excitonic enhancement is ac- 
companied by a region of negative differential conduc- 
tance. 

Color-scale plots of the current as a function of bias 
voltage V = {hl — Mfl.)/ e an d gate voltage V g = —e d /e 
are shown in Fig. 2J The excitonic interaction leads to 
a polaron shift, A, which renormalizes the dot charing 
energy, U — >• U = U — 2A, and thus reduces the size of 
the Coulomb diamonds. 

As discussed in Sec. |VJ asymmetric excitonic dot-lead 
couplings change the behavior. In Fig. [5] we demon- 
strate this behavior, presenting results calculated for pa- 
rameters appropriate to a nanotube (M = 4 channels 
and K c ~ 0.7) but with an excitonic interaction only 
to the left lead, so that Zl ~ 1 and Zr ~ 0, imply- 
ing Yl ~ 0.75 and Yr ~ 1.5. In this case tunneling to 
the left lead is enhanced and tunneling to the right lead 
is suppressed. To understand the consequences for the 
I-V characteristics we refer to Eq. ([55)) for the steady- 
state current. Near the lowest threshold we may ne- 
glect the 1— >2 transitions and simplify the expression to 
(I L ) = 2e(ft£_> 1 ftf_ +0 -ftf_> ^i)/(2ft^i+rci-fo). 
We see that when the bias voltage reaches the positive 
threshold eV = +2£d, only the first term contributes be- 
cause TZn-i ~ in this case. Since the rate TZP.n is a 
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(a) 

I L (eh L ) 




FIG. 5: (a) Current-voltage characteristics and (b) transition 
rates for asymmetric excitonic dot-lead couplings, Yl = 0.75 
and Yr = 1.5. We assume T = id/WO, v c A — 15id and as 
in Fig. [2] choose U = led- Voltages are given in units of the 
onsite energy. 



smooth function at eV = +2e ( i and the denominator has 
no singularities, the current is essentially proportional to 
eT^o-n- However, when the bias voltage reaches the neg- 
ative threshold eV = —2id, only the second term con- 



tributes because 72-o->i — hi that case. Accordingly 
the current is asymmetric and shows step-like features at 
positive bias voltages but a suppressed onset at negative 
bias voltages. In the asymmetric case the "Luttinger" 
renormalization of the tunneling amplitude is different 
for tunneling into the two leads [cf. Eq. (|60])]. also con- 
tributing to the asymmetry. 



VII. CONCLUSIONS 

In summary, we have investigated the effect of a 
Coulombic dot-lead coupling on transport through a 
quantum dot coupled to two semi-infinite Luttinger- 
liquid leads. The electronic tunneling has been described 
within a master-equation approach that treats the dot- 
lead hybridization to lowest nonvanishing order. We have 
found that the excitonic dot-lead interaction may en- 
hance transport in the vicinity of the Coulomb-blockade 
threshold, which is in contrast to the power-law suppres- 
sion of the electronic tunneling if this interaction is not 
included. However, the effects are in general less pro- 
nounced than for the side-coupled situation considered 
in our previous work,— because in the present two-lead 
situation both the excitonic and orthogonality effects arc 
present. 

The paper raises several questions for future research. 
A treatment of the electronic tunneling to all orders in 
the hybridization would be desirable. Moreover, it would 
be interesting to study the effect of the Coulombic dot- 
lead coupling on transport in the Kondo regime, where 
electronic tunneling is dominated by dot-lead exchange 
scattering processes. Finally, the consequences of an in- 
teraction between the two leads should be explored. 
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